Numerical approach to a model for quasistatic damage with spatial BV-regularization

We address a model for rate-independent, partial, isotropic damage in quasistatic small strain linear elasticity, featuring a damage variable with spatial BV-regularization. Discrete solutions are obtained using an alternate time-discrete scheme and the Variable-ADMM algorithm to solve the constrained nonsmooth optimization problem that determines the damage variable at each time step. We prove convergence of the method and show that discrete solutions approximate a semistable energetic solution of the rate-independent system. Moreover, we present our numerical results for two benchmark problems

Fully discrete approximation of rate-independent damage models with gradient regularization

This work provides a convergence analysis of a time-discrete scheme coupled with a finite-element approximation in space for a model for partial, rate-independent damage featuring a gradient regularization as well as a non-smooth constraint to account for the unidirectionality of the damage evolution. The numerical algorithm to solve the coupled problem of quasistatic small strain linear elasticity with rate-independent gradient damage is based on a Variable ADMM-method to approximate the nonsmooth contribution. Space-discretization is based on P1 finite elements and the algorithm directly couples the time-step size with the spatial grid size h. For a wide class of gradient regularizations, which allows both for Sobolev functions of integrability exponent r ∈ (1, ∞) and for BV-functions, it is shown that solutions obtained with the algorithm approximate as h → 0 a semistable energetic solution of the original problem. The latter is characterized by a minimality property for the displacements, a semistability inequality for the damage variable and an energy dissipation estimate. Numerical benchmark experiments confirm the stability of the method

Fully discrete approximation of rate-independent damage models with gradient regularization

This work provides a convergence analysis of a time-discrete scheme coupled with a finite-element approximation in space for a model for partial, rate-independent damage featuring a gradient regularization as well as a non-smooth constraint to account for the unidirectionality of the damage evolution. The numerical algorithm to solve the coupled problem of quasistatic small strain linear elasticity with rate-independent gradient damage is based on a Variable ADMM-method to approximate the nonsmooth contribution. Space-discretization is based on P1 finite elements and the algorithm directly couples the time-step size with the spatial grid size h. For a wide class of gradient regularizations, which allows both for Sobolev functions of integrability exponent r ∈ (1, ∞) and for BV-functions, it is shown that solutions obtained with the algorithm approximate as h → 0 a semistable energetic solution of the original problem. The latter is characterized by a minimality property for the displacements, a semistability inequality for the damage variable and an energy dissipation estimate. Numerical benchmark experiments confirm the stability of the method

Approximation schemes for materials with discontinuities

Damage and fracture phenomena are related to the evolution of discontinuities both in space and in time. This contribution deals with methods from mathematical and numerical analysis to handle these: Suitable mathematical formulations and time-discrete schemes for problems with discontinuities in time are presented. For the treatment of problems with discontinuities in space, the focus lies on FE-methods for minimization problems in the space of functions of bounded variation. The developed methods are used to introduce fully discrete schemes for a rate-independent damage model and for the viscous approximation of a model for dynamic phase-field fracture. Convergence of the schemes is discussed

Approximation schemes for materials with discontinuities

Damage and fracture phenomena are related to the evolution of discontinuities both in space and in time. This contribution deals with methods from mathematical and numerical analysis to handle these: Suitable mathematical formulations and time-discrete schemes for problems with discontinuities in time are presented. For the treatment of problems with discontinuities in space, the focus lies on FE-methods for minimization problems in the space of functions of bounded variation. The developed methods are used to introduce fully discrete schemes for a rate-independent damage model and for the viscous approximation of a model for dynamic phase-field fracture. Convergence of the schemes is discussed

Primal-dual gap estimators for

The primal-dual gap is a natural upper bound for the energy error and, for uniformly convex minimization problems, also for the error in the energy norm. This feature can be used to construct reliable primal-dual gap error estimators for which the constant in the reliability estimate equals one for the energy error and equals the uniform convexity constant for the error in the energy norm. In particular, it defines a reliable upper bound for any functions that are feasible for the primal and the associated dual problem. The abstract a posteriori error estimate based on the primal-dual gap is provided in this article, and the abstract theory is applied to the nonlinear Laplace problem and the Rudin–Osher–Fatemi image denoising problem. The discretization of the primal and dual problems with conforming, low-order finite element spaces is addressed. The primal-dual gap error estimator is used to define an adaptive finite element scheme and numerical experiments are presented, which illustrate the accurate, local mesh refinement in a neighborhood of the singularities, the reliability of the primal-dual gap error estimator and the moderate overestimation of the error

Approximation schemes for materials with discontinuities

Damage and fracture phenomena are related to the evolution of discontinuities both in space and in time. This contribution deals with methods from mathematical and numerical analysis to handle these: Suitable mathematical formulations and time-discrete schemes for problems with discontinuities in time are presented. For the treatment of problems with discontinuities in space, the focus lies on FE-methods for minimization problems in the space of functions of bounded variation. The developed methods are used to introduce fully discrete schemes for a rate-independent damage model and for the viscous approximation of a model for dynamic phase-field fracture. Convergence of the schemes is discussed